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ABSTRACT 

O ■ 

\Q • The two particle inclusive cross section for the reaction p + p — > 

7 + c + X is studied in perturbative quantum chromodynamics at 
^ ■ order 0{a 2 s ), for large values of the transverse momentum of the 

Qh! prompt photon and charm quark. Two different techniques are used 

o ■ in performing the phase-space integrals; the first is purely analyti- 

D \ cal, and the second is a combination of analytic and Monte Carlo 

& • integration methods. The second, more versatile technique facili- 

tates imposition of photon isolation restrictions and other selections 
■ of relevance in experiments. Differential distributions are provided 

rS | for various observables, and a comparison is made with preliminary 

data from the CDF collaboration. 
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1. Introduction 

The observation of a photon earring a large transverse momentum, pt, among the final 
state products of a high energy hadronic reaction has long been regarded as a powerful tool 
for gathering information about the short distance dynamics of these reactions, primarily 
because photons couple directly to quarks via the point-like electromagnetic interaction. 
Single prompt photon production has been studied extensively, both theoretically and ex- 
perimentally, and has provided useful information about hadronic structure, in particular 
the gluon structure function of the proton. Data are beginning to appear from the study of 
the production of a prompt photon (7) in association with a heavy quark whose transverse 
momentum balances a substantial portion of that of the photon. This two-particle inclusive 
reaction is particularly interesting because it offers the possibility of a detailed study of the 
underlying QCD dynamics such as, e.g., rapidity correlations. In addition, it may provide a 
direct measurement of the charm content of the proton. 

We report here on two next-to-leading order perturbative QCD calculations we have done 
of the reaction p + p^7 + c + Xat high energy 1,2 \ In these calculations two different 
techniques are used in performing the phase-space integrals. In the first, purely analytical 
techniques are used. In the second approach, we use a combination of analytical and Monte 
Carlo techniques, which is more flexible and allows implementation of isolation cuts and 
other experimentally relevant selections. To warrant use of perturbation theory and the 
massless approximation, we limit our considerations to values of transverse momenta of the 
photon and charm quark p^ c > 10 GeV. 

Only one direct hard scattering subprocess contributes in leading order: the quark-gluon 
Compton subprocess gc — > 7c. The initial charm quark and the initial gluon are constituents 
of the initial hadrons. In addition, there is a leading order fragmentation process in which the 
photon is produced from quark or gluon fragmentation, e.g., gg — > cc followed by c — > jX, 
or qc — > qc followed by q — > 7. At next-to-leading order in QCD, several subprocesses 
contribute to the 7 + c + X final state: gc — > gey, gg — > eery, qq — > ccy, qc — * qcy, qc — > qcy, 
cc — > eery, and cc — > CC7. A full next-to- leading order calculation requires the computation 
of the hard-scattering matrix elements for these two-to-three particle production processes 
as well as the one- loop 0(a 2 ) corrections to the lowest order subprocess gc — > 7c. At 0(a 2 s ) 
there are, in addition, fragmentation processes in which the hard-scattering two-particle 
final-state subprocesses c + g ^ 7 + c, c + c 7 + g and q + q — > 7 + g are followed by 
fragmentation processes c — > cX, in the case of the first subprocess, and g — > cX in the 
cases of the last two. These should be included because we factor the collinear singularities 
in the corresponding three-body final-state processes into non-perturbative fragmentation 
functions for production of a charm quark from a particular parton. 
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We are interested in the two-particle inclusive differential cross section, E 7 E c da / d 3 p y d 3 p c , 
where (E, p) represents the four- vector momentum of the 7 or c quark. For each contributing 
subprocess, this calculation requires integration over the momentum of the unobserved final 
parton in the two-to-three particle subprocesses (g, c,q, or q). Collinear singularities must 
be handled analytically by dimensional regularization and absorbed into parton momentum 
densities or fragmentation functions. 

2. Analytic Calculation 

In our purely analytical calculation 1 , we must integrate over enough of the phase-space 
analytically in order to cancel all soft and collinear singularities. This procedure imposes 
technical limitations such that it is not possible to obtain the fully differential two-particle 
cross section. We calculate instead the cross section da / dpTpdy 1 dz , where the variable z is 
defined by 

This definition indicates that z contains information on the transverse momentum of the 
charm quark, but no information on its rapidity y c . Moreover, it is not possible to implement 
isolation cuts on the photon in the purely analytic calculation. Nevertheless the calculation 
is valuable for obtaining the relative importance of the various subprocess contributions and 
the overall the size of the cross section, as well as for checking aspects of the more versatile 
Monte Carlo calculation. 

We calculate the cross section da 7 'dp^dxf 'dz . Whenever the charm quark and photon 
have balancing values of pt then z — 1. This occurs 

1. for the leading order process eg — > 7c, 

2. for the virtual gluon exchange processes, 

3. whenever a gluon becomes soft in the three-body processes, or 

4. when a third parton becomes collinear to the beam in a three-body process. 

The point z — 1 is associated with various soft and collinear poles, and in order to expose 
these poles all phase-space integrals are carried out in 4 — 2e dimensions. The poles are 
exposed by expanding the results in (1 — z). The analytical cross section is expressed in terms 
of distributions in (1 — z) such as 1/(1 — z)+ and 6(1 — z). In addition whenever the final state 
photon becomes collinear to a final state quark, the resulting singularity must be exposed in 
a similar way before it is absorbed into the photon fragmentation function. This singularity 
occurs at another value of z, (z = z-y). This means that there are also distributions in (z — Zi). 
The cross section is singular at z — 1 and z — z±, but the singularities are integrable through 
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Fig. 1. Cross section as a function of z at y 1 = in (a) leading order and (b) next-to-leading showing 
the dominant contributions. 



the use of 'plus'-distributions. The cross section is presented in terms of integrals over finite 
regions of z such as 



da 1 r z+ ^r da , 

~ Az dpldy^dz' Z ' 

dz'. (2) 



d-p^dxpdz Az Jz~^ d-p^dxpdz 
da r z b da 



dpTpdy^ Jz a dpJpdy^dz' 

Our results are presented at a center-of-mass energy y/s = 1.8 TeV. The analytical 
results are obtained with the GRV parametrization of the proton parton densities, and all 
renormalization/factorization scales are taken as /i — pj. Up to 30% differences are observed 
if the CTEQ3M parton distributions are used. We sum over charm and anticharm production 
throughout. In Fig. la the net lowest order contribution is shown as a function of z at p7p = 15 
GeV and for a bin size Az = 0.2. The lowest order cross section is made up of the lowest 
order direct term eg — * 7c, which is proportional to 5(1 — z) and provides the peak at z — 1, 
plus the various photon fragmentation contributions which contribute in the region z > 1. 
The most significant feature of this curve is that there is no contribution to the cross section 
in the region z < 1. This unrealistic prediction shows the inadequacy of the lowest order 
predictions. 

Figure lb shows the distribution in z predicted by the next-to-leading order calculation. 
The next-to-leading order contributions serve to lower the peak at z — 1, and they broaden 
the distribution. The cross section is finite at all values of z, closer to the situation observed 
in experiments. In addition, in Fig. lb we display contributions from the most important 
subprocesses. The eg initiated process dominates the cross section, but there are important 
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contributions from the gg and cq initiated process in the low pj. region. 

3. Monte Carlo Calculation 

The combination of analytic and Monte Carlo techniques used here to perform the phase 
space integrals is documented and described elsewhere 2 . The analytic/Monte Carlo method 
allows for the calculation of many different observables without much extra effort and for 
the imposition of experimental cuts. In collider experiments a photon is observed and its 
momentum is well measured only when the photon is isolated from neighboring hadrons. In 
our calculation, we impose isolation in terms of the cone variable R: 

y/(Ay)* + (A0)2 < R. (3) 

In Eq. (|3p, Ay (A0) is the difference between the rapidity (azimuthal angle in the transverse 
plane) of the photon and that of any parton in the final state. The photon is said to be 
isolated in a cone of size R if the ratio of the hadronic energy in the cone and the transverse 
momentum of the photon does not exceed e = 2GeV/pp. We show distributions for the 
choice R = 0.7 typical of current experiments. 

The structure of the QCD hard-scattering matrix element produces positive correlations 
in rapidity at collider energies. To examine these correlations more precisely, in Fig. 2a, 
we display the differential cross section in y c , for two intervals of y 1 in the forward rapidity 
region. These distributions show that the typical rapidity of the charm quark follows that 
of the photon, and thereby confirms the positive correlation between the rapidities of the 
photon and charm quark. In Fig. 2b we compare our results to CDF data 3 for photon plus 
[i^ production. The three upper points are obtained 4 from the Monte Carlo event generator 
Pythia whereas the lower point is that given by our theoretical calculation. The Pythia 
results lie substantially below the data whereas our results are higher than the data but 
somewhat closer to it. 

4. Conclusions 

We presented the results of two calculations of the inclusive production of a prompt 
photon in association with a heavy quark at large values of transverse momentum. Both 
analyses are done at next-to-leading order in perturbative QCD. Our results agree quanti- 
tatively as they should, but the combination of analytic and Monte Carlo methods is more 
versatile. We provide differential cross sections in transverse momenta and rapidity, includ- 
ing photon isolation restrictions, that should facilitate contact with experimental results at 
hadron collider energies. We show that the study of two-particle inclusive distributions, with 
specification of the momentum variables of both the final prompt photon and the final heavy 
quark, tests correlations inherent in the QCD matrix elements and should provide a means 
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Fig. 2. (a)Cross section da/dy c as a function of the rapidity of the charmquark forp+p — > 7 + c+X 
at y/s — 1.8 TeV. The solid curve shows the result when the photon rapidity is restricted to 1.0 < 
y 1 < 2.0, and the dashed curve displays the result for 2.0 < y 1 < 3.0. (b) Ratio of the measured 
cross section to that predicted by theory for various final state charm decay products. 

for measuring the charm quark density in the nucleon. A comparison of our next-to-leading 
predictions with the preliminary CDF data shows reasonable agreement. 
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